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Abstract. The null geodesies of a Schwarzschild black hole are studied from a 
dynamical systems perspective. Written in terms of Kerr-Schild coordinates, the 
null geodesic equation takes on the simple form of a particle moving under the 
influence of a Newtonian central force with an inverse-cubic potential. We apply 
a McGehee transformation to these equations, which clearly elucidates the full 
phase space of solutions. All the null geodesies belong to one of four families of 
invariant manifolds and their limiting cases, further characterized by the angular 
momentum L of the orbit; for \L\ > \Lc\, (1) the set that flow outward from the 
white hole, turn around, then fall into the black hole, (2) the set that fall inward 
from past null infinity, turn around outside the black hole to continue to future 
null infinity, and for \L\ < \Lc\, (3) the set that flow outward from the white hole 
and continue to future null infinity, (4) the set that flow inward from past null 
infinity and into the black hole. The critical angular momentum Lc corresponds 
to the unstable circular orbit at r = 3M, and the homoclinic orbits associated 
with it. There are two additional critical points of the flow at the singularity 
at r = 0. Though the solutions of geodesic motion and Hamiltonian flow we 
describe here are well known, what we believe is a novel aspect of this work is the 
mapping between the two equivalent descriptions, and the different insights each 
approach can give to the problem. For example, the McGehee picture points to a 
particularly interesting limiting case of the class (1) that move from the white to 
black hole: in the L — > oo limit, as described in Schwarzschild coordinates, these 
geodesies begin at r = 0, flow along t = constant lines, turn around at r = 2M, 
then continue to r = 0. During this motion they circle in azimuth exactly once, 
and complete the journey in zero affine time. 
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1. Introduction 

The Schwarzschild metric, describing a static, non-rotating black hole solution to the 
Einstein field equations, was discovered within a year of Einstein's publication of the 
theory of general relativity. Though it had to wait till the 1960's before the full nature 
of the metric was truly uncovered, it has nevertheless been well studied for nearly a 
century. One of the more important tools in this regard is understanding the geodesic 
structure of the spacetime; in particular, timelike and null geodesies characterize the 
paths that a freely moving test particle P can follow within the spacetime. The 
geodesic equations of motion for P are a set of second order, ordinary differential 
equations, describing the evolution of the coordinates of P as a function of an affine 
parameter. In Schwarzschild these equations are integrable, and the solutions have 
been known for a long time (see any standard text on general relativity). Nevertheless, 
the solutions exhibit a sufficiently rich set of dynamics that new insights into them are 
still being garnered. For example, in the past couple of decades dynamical systems 
methods have been brought to bare on the geodesic equations, and perturbations 
thereof. A couple of the interesting results have been a new taxonomy of orbits based 
on the subset of orbits that are periodic, emphasizing the importance of the homoclinic 
orbits that asymptotically approach the unstable branch of circular orbits [H [H [3] , 
and that a generic perturbation of the geodesic flow possesses a chaotic invariant 
set [l[51ll[71IHlini[TO|. 

In this paper we describe a new approach to understand the phase space 
of geodesic orbits, by applying a transformation due to McGehee [11]. This 
transformation is designed to resolve the singularities that formally appear in the 
Newtonian equations of motion when particles interacting through a central force 
collide. The transformation "blows up" collision into an invariant manifold for the 
flow of the transformed differential equations. This allows the motion of P to be 
studied near collision, uncovering interesting dynamics. To apply these methods, 
we map the geodesic equation, written in Cartesian-like Kerr-Schild 12 (or ingoing 
Eddington-Finkelstein [l3j[T4j) coordinates, to a central force problem. This is trivial 
for the case where P has zero mass (null geodesic), which we focus on in this paper, 
though the methods can be generalized. 

The rest of the paper is laid out as follows. In Section [2] we introduce the 
Schwarzschild metric in the standard Schwarzschild coordinates, give an overview of its 
Penrose diagram, and qualitatively discuss the null geodesies of the geometry. This will 
provide a point of reference as we map the geometric picture to a dynamical systems 
description within McGehee coordinates. In Section [3] the differential equations for 
the motion of P are given, in Schwarzschild and Kerr-Schild coordinates. In Section!?] 
we briefly summarize the results of the McGehee method for understanding the flow 
of a general class of Hamiltonian systems corresponding to central force problems. In 
Section [5] we restrict to the particular system that maps to the null geodesic structure 
of Schwarzschild. We show that the phase space flow can be subdivided into four 
families of invariant manifolds and their limiting cases. There are four critical points 
of the flow — the black hole and white hole singularities at r = 0, and the two (with 
angular momentum L = ±Lc) unstable circular orbits at r = 3Af . Though much 
of the understanding of geodesic motion and the classical Hamiltonian system are 
individually well known, we believe the novel aspect of this work is the (sometimes 
non-trivial) mapping between the two equivalent descriptions, and the insight one gives 
to the other. One particularly interesting example is the limiting case of geodesies that 
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in the dynamical systems description flow "directly" from the the white to black hole: 
the collective set of these geodesies trace out the interior region of the white/black hole, 
flow along Schwarzschild t = const, lines, and complete exactly one orbit A$ = 27r 
along the journey from white hole to black hole singularity. Also, the standard affine 
parameter integrates to zero along these curves. To obtain a finite afBne length requires 
rescaling it by an angular momentum dependent quantity, which diverges in the limit. 
Details of all this are discussed in Section [5l Finally, we conclude in Section [6] with a 
summary, and discussion of possible future extensions and applications. Throughout 
we use geometric units, where the speed of light c = 1 and Newton's constant G — 1. 

2. The Schwarzschild Geometry 

The Schwarzschild metric, describing a non-rotating black hole of mass M , has the 
following line element in standard (spherical polar) Schwarzschild coordinates: 

ds"^ = gai3dx"dx^ = -(1 - 2M/r)dt^ + (1 - ^M/ry^dr"^ + r'^{de'^ + sin^ Odcf)'^), (1) 

where gap is the metric tensor, and we use the Einstein summation convention 
where repeat indices are summed over. Figure [1] shows a Penrose diagram of the 
Schwarzschild spacetime. For readers not familiar with this diagram, we will briefly 
review its salient features. The Penrose diagram is a conformal compactification of 
the maximal analytic extension of the Schwarzschild metric, designed to highlight 
the causal structure of the spacetime. Here we project out the {9, (p) coordinates, so 
a point on the diagram represents a two-sphere of area Aur^ . Radial (zero angular 
momentum) null geodesies are straight lines angled at ±45° relative to the horizontal. 
Any causal curve that could be associated with a particle trajectory (whether timelike, 
null, geodesic or not) has a slope > 45° everywhere along its projected curve on the 
diagram. There are two identical, causally disconnected, asymptotically flat regions 
reached as r — >■ cxd; without loss of generality we focus on the right hand region on 
the diagram. Similarly, there are two singularities in the geometry at r = 0: the white 
hole singularity occurring to the past of any event in the spacetime, and the hlack hole 
singularity that is within the causal future of any event. 

A null curved can be parameterized by an affine time cr, unique up to a constant 
scaling and translation. Null curves that originate at r = cx) (with a — > —oo) come 
from a region called past null infinity 3~ , and those that return to r = cx) (as a -\-oo) 
end at future null infinity 5+ . The event horizon of the black hole is formally defined 
as the boundary of the causal past of 5"*", which on the Penrose diagram is the line 
labeled r — 2M,t = oo. Null curves that come from the white hole "begin" at r = 
with finite a, and those that cross the event horizon reach r = in finite affine time. 
Note that the Schwarzschild time coordinate t is singular at r = 2M. Further, note 
that r and t switch "character" at r = 2M: outside, r is spacelike, t timelike, and 
vice-versa inside. In the black hole, time as measured by r flows to smaller r (the 
opposite in the white hole). Thus, the inevitability of encountering the singular at 
r = for any observer crossing the horizon is evident in the diagram: r = is no 
longer a "place" that can be avoided, rather it is a "time" that will happen for any 
such observer. Finally, note that the compactification severely distorts the spacetime 
at the corner points of the diagram — that r = 0, r = 2M and r — oo touch in the 
diagram is purely an artifact of the compactification. 




Figure 1. A schematic Penrose diagram of the Schwarzschild solution, with 
representative null geodesies from each of the invariant manifolds; see Figure [2] 
for the corresponding diagram in McGehee coordinates. 



3. Geodesic Equations in the Schwarzschild Spacetime 

We begin by describing the geodesic equations for causal particles; later we restrict to 
the null case. The geodesic equation for a particle P in parametric form x°'{a) is 

+ T^^x^x'' = (2) 

where the over-dot ' = d/ da and V^^ is the metric connection. This is formally a set 
of second order ordinary differential equations, though due to the symmetries of the 
Schwarzschild metric, and the normalization condition that gai3X°'x^ equals (— 1) 
for null (timelike) geodesies, one obtains the following first integrals of motion of ([2]) 
in Schwarzschild coordinates ([Ij: 

i^E{l- 2M/r)-i 

^ ^2^4 _ (-^2 _ 2Mr){fi^r^ + K) 
r'^Q^ = K -L^ sin"2 q 

,P ^ L{r^ sin^ 0)-\ (3) 

where /i^ is the rest mass of P, E its energy, L its angular momentum about the axis 
sin0 = 0, and K is Carter's constant of motion. Note that the second equation above 
for f can be written as a Newtonian central force problem 

£; = ir^ + l/(r), (4) 



with equation of motion 



f = —dV/dr, 



(5) 
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where E = & /2 and the effective potential V{r) is given by 

Kerr-Sehild coordinates 

The Schwarzschild metric in Cartesian Kerr-Schild coordinates (A,x,y, z) is 

ds^ = -d\^ + dx^ + dy^ + rfz^ + 2Mr^^{xdx + ydy + zdz + rdXf , (7) 

related to Schwarzschild coordinates by the set of transformations 

a; = r sin cos 0, y = r sin sin 0, z~rcos9 (8) 

and 

dX = dt^ dr. (9) 

One advantage of Kerr-Schild coordinates over Schwarzschild coordinates, as is evident 
from (O, is that they are regular across the event horizon at r = 

The geodesic equations in Kerr-Schild coordinates can be written as [l5] 

X = 2MKr-^{E-r)-^ +E (10) 

ifc = -SMKxkr^^ - fi^Mxkr'^, (11) 

where for the sake of notation, k — 1,2,3, with xi — x,X2 ~ y,xz — z. Due to the 
spherical symmetry of the Schwarzschild spacetime, without loss of generality we can 
restrict attention to planar motion. For simplicity we choose z = 0, corresponding to 
6 — 7r/2, for which K ^ In the expressions below we will replace K with L^, 

though one can consider them to be valid for motion in any plane if L is re-interpreted 
as the angular momentum relative to an axis orthogonal to the orbital plane. We will 
also now limit the discussion to null geodesies, for which — 0, and only focus on 
the coordinate flow pT|) of the geodesies: 

X = -SML'^xr-^, (12) 



where x — {xi,X2) G R^, and r = |x| = y/x^ + x^. 

The corresponding Hamiltonian for this system of equations is 

i?-i/(x,x)-i|x|2-L2M|x|-3, (13) 

where E — E'^ /2 is a positive constant of motion for any particle trajectory. 

f Technically, ^ is only regular across the black hole event horizon, and not the Cauchy horizon at 
r = 2M of the white hole. The reason is A has been chosen so that along an r = const, surface, A 
coincides with the retarded time v of an ingoing radial null curve coming from S~, and v — > — oo 
as the white hole horizon is approached on the Penrose diagram. These coordinates are thus also 
referred to as ingoing Eddington-Finkelstein coordinates. If instead we had chosen a coordinate 
X = dt — 2M/{r — 2M)dr, then A would coincide with advanced time u of an outgoing radial photon 
along an r = const, surface (outgoing Eddington-Finkelstein coordinates), and the metric would be 
regular at the white hole horizon, but singular at the black hole horizon for similar reasons. However, 
at the end of the day one arrives at the same geodesic equation of interest II12I I for the spatial 
Cartesian coordinates as a function of affine time [a;(cr), y{cr), z(a-)], which is well defined across both 
horizons. 
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It is noted that for each value /i > of the energy the motion of P of mass zero 
hes on the three-dimensional energy surface 

S = {(x,x) e i?^|i/(x) = > 0}. (14) 



To facilitate direct application of prior work on Newtonian central force problems 
of the form (|12p. it is useful to perform the following dimensionful rescaling of the 
afhne parameter ct to ^ 

a-^^. (15) 

This transforms ([T^ to 

^ = -3xr-, (16) 
In the following section we consider the solutions x(^) ~ {xi{S^), X2{£,)) of (fTB)) . 

4. Transformation to McGehee Coordinates and Blowup of Collision 

We consider general central force fields with potential 

C/(x) = -|xr", (17) 

a > 0, later restricting to the case a = 3 of ^TE\i . The system of differential equation 
describing the motion of a single particle with this potential is given by 

x = -9xC/(x) = -a|x|-"~2^. (18) 

It is convenient to write this as the first order system 
X = y, 

y = -a|xr"-2^. (19) 
This is a Hamiltonian system with Hamiltonian function 

77(x,y) = i|yp-|xr", (20) 

which is the total energy of the particle and is conserved along solutions of p^ . where 
(x,y)ei?4. 

The general description of the orbit structure for (jl9p was carried out in |11| . 
with a special emphasis on the motion of the particle near collision. The general 
approach taken is to find a change of coordinates which have the effect of blowing 
up the collision, corresponding to x = 0, into an invariant manifold with its own 
flow. When this is done, the dynamics of the particle near collision can be completely 
understood and solutions tending towards collision are asymptotic to this manifold. 
This change of coordinates also gives the global flow for the differential equations. 

Set X — (x, y), and consider a solution X(^) ^ (x(^),y(f)) for (fTOl) with an 
initial condition X(0). The standard existence and uniqueness theorems of differential 
equations guarantee that X(^) can be uniquely determined and defined over a maximal 
interval (^~,C^), where — oo < ^~ < < ^+ < +00. 
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Definition If ^+ < oo, then X(^) ends in a singularity at If > —oo then 
X(^) begins in a singularity at In either case, ^* = ^+ or is said to be a 
singularity of the solution X(^). 

The following result is proven in [IJJ: 

Let X(^) be a solution of ([T9l) with a singularity at Then this singularity is due 
to collision. That is, x(^) — > as z — 

It follows from (PO)) that for a collision solution, |y| — oo as 2; — Hence, 
for (|19p the only solutions that are singular are collision solutions, either ending or 
beginning in collision. There are several methods available to study collision, the 
dynamics of solutions near it and to understand if a collision solution can be extended 
through the collision state in a smooth fashion. Here we present a brief summary of 
blowing up collision for (|19l) for arbitrary a > 0, and then apply that to understand 
the global flow for the case of a = 3. The details are in 

We set (3 — ^, 7 = [1 + /3)^^ for sake of notation. Also, it is convenient to use 
complex coordinates and identify the real plane with the complex plane . Then, 
we can consider x to be a vector in or a complex number in C^. The McGehee 
coordinates are given by the transformation T of (xi, X2, yi, j/2) to (r > 0, 9, w, v), 

x = Pe'^ (21) 
y = f-l^^ [v + iw)e'-^^ 
and a transformation of the afhne variable f , 

d£, = fdr. (22) 

The system (fT9|) is transformed into 
f' = + l)fv, 

9' = w, (23) 

w' — {l3 — l)wv, 

v' = +P{v^ - 2), 

where prime denotes differentiation wrt r. In complex notation, the angular 



momentum for (|19p is given by 

l](x,y) =3(xy). (24) 

We fix the energy H to the constant value h and O to the constant value c (not to be 
confused with the speed of light). T transforms H — h and Vl — c into 

+v'^ -2^ 2hr°''' (25) 
= c. (26) 
We define the constant energy manifold, 

M{h) = {(f , 9, w, v) e R^\f>0,H = h}. (27) 

We define the collision set corresponding to collisions for System as 



N = {{f;0,w,v) e M(/i)|f = 0}. 



(28) 
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On account of (pSj) . (PB)) and ([25)1 can be written as 

N = e M(/i)|f = 0,t;2 2,u; = 0} /3 > 1, (29) 

N = G M(/i)|f = 0,t;2 + = 2,w = c} ^ = 1, (30) 

N = {{f,e,w,v) e m{h)\f = 0,w^ = 2,c = 0} < ^ < 1. (31) 

As is proven in [TT] , 



N is an invariant manifold for the vector field defined by System ([23|) . Collision orbits 
approach N asymptotically as r — > ±cx). 

Definition N is called a blow up of the collision f = on M(/i). 

It is also shown in [TI] that System reduces to the system consisting of the last 
two equations of 

= (/3 - l)u;w, (32) 
v' = _^^(t;2 - 2). 
This system has the integral 

A(w,z;) = |w;|"|v2+w2_2|i-'^. (33) 

In summary, to understand the flow of (IT^ on M(/i), we can use System (P^ . where 
the collision manifold is given by N. 

5. McGehee Flow of Schwarzschild Null Geodesies 

We now consider the case of interest given by the Schwarzschild problem that we 
showed reduced to System p^ . This corresponds to a = 3, /3 — 3/2, 7 — 2/5. The 
flow for System (|32|) is depicted in Figure [2l Incidentally, this figure captures the 
qualitative behavior for all values of /3 > 1. It turns out that the fiows for /3 < 1 
are drastically different, in particular the Kepler problem corresponds to /3 = 1/2. 
Over the next several paragraphs we will describe the fiow in some detail, though 
first it would be helpful to identify a few relations between the coordinates {w, v) and 
physical coordinates, as well as establish the map between the constants of motion 
{h,c) of solutions to System p9|) and the physically relevant constants L,E and M. 

Relationships between constants and coordinates 

Using (PT|) and (^51) , and the angular momentum of the geodesic in the physical picture 
L — xy — yi, one finds that the angular momentum c of the System ([T^ evaluates to 

c = (34) 

This rather unusual relationship is due to the scaling P3|) between the time parameters 
of the two descriptions. Likewise, it is straightforward to show that the energy h of 
the System ([TO]) relates to physical constants via 

From (|2ip . the relation between areal radius r and f for a = 3 is 

r = (36) 
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Figure 2. The Schwarzschild null geodesic flow in McGehee Coordinates; see 
Figure [T] for the corresponding Penrose diagram. The point P~ {P^) here 
corresponds to the black hole (white hole) singularity at r = 0. For clarity we 
only show the half plane w > 0, corresponding to positive angular momentum L 
orbits; the flow in u) < for negative angular momentum orbits looks identical to 
that of «) > 0, but reflected about the v axis. 



and this together with and (|26p gives 

/ r 

(37) 




where the plus (minus) sign corresponds to positive (negative) physical angular 
momentum L. Thus vertical lines of constant w on the McGehee diagram correspond 
to surfaces of constant areal radius r. Putting all the above together, with ([25l) . one 
finds that all trajectories projected onto the {w,v) plane are characterized by the 
following polynomial with one free parameter Q > 0: 

= 2 - + Qw^, (38) 

Q ^ -j^- (39) 

Finally, note from ([25]) that v ~ implies r' ~ 0, which translates to f = 0; thus a 
trajectory that crosses v = corresponds to a geodesic that has a turning point in its 
radial motion. 

Particle Flow 

Figure [5] is a projection of the full flow in {f,d,w,v) coordinates to the {w,v) plane. 
Though as seen in (p3)) . the two differential equations for {w',v') do not depend on 
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(f, 6) and can be solved separately as an independent system. Once we know [w, v) 
then the coordinates {f,0) can be easily computed for each {w,v) by (1^51) . {6 and 
r) represent polar-like coordinates for P, while {w, v) can be viewed as velocity-like 
coordinates. System implies that 9{t) either increases, w > 0, or decreases 
w < 0. This is just a cycling motion about the origin f = in position coordinates 
f, 0. While this cycling motion is occurring, r increases, u > 0, or decreases, v < 0. 
This is analogous to the projection on the Penrose diagram of the flow in (r, t, 6, (f)) to 
the (conformally compactified) (r, t) plane, and cyclic motion for K = geodesies 
corresponds to increasing (/>, L > 0, or decreasing 0, L < 0. 

The flow curves {'w{t),v(t)) in Figure [2] can be viewed as invariant manifolds. 
They foliate the {w, w)— plane. From (l23l) one can see that the flow has critical points 

at {w,v) = (±-\/3, 0), where r — a/3, and 9 = ±-\/3t (to within a constant phase) 
respectively. These correspond to the unstable (hyperbolic) circular periodic orbits at 
r = 3M. The periodic orbits exist for each /i > on each energy surface M(/i) . 

Excluded region 

The flow on M(/i) for ([32|) projects into the set {w'^ + v'^ < 2} for each /i < 0, to 
{w^ + — 2} for h = 0, and to {w^ + > 2} for h > 2. Hence, since P moves on 
E, defined by (fT4| . where h > 0, then we need not consider the interior of the disc, 
{w^ + < 2}. So, in Figure [5] the region of interest are all points {w^ + > 2}. 

The collision set 

The collision set N for a = 3 reduces to two critical points of the flow, p+ ~ (0, -\/2) 
and p~ = (0,— a/2)- These are unstable hyperbolic points. The flow tends to p^ as 
T — > — oo and to as r — > +oo. In the full phase space, these critical points have a 
constant value of 9, and f = 0. On the Penrose diagram, p~ {p'^) corresponds to the 
black hole (white hole) singularity at r = 0. 

Invariant manifolds of the flow 

There are invariant manifolds connecting and p^ . Two flow from p^ to and 
two flow from to p^ . These we label W^^^ and , respectively; trajectories 
within this flow have physical angular momentum-squared = 27 E^M^. As P moves 
along, e.g. , it can be viewed in position space as cycling about the origin while 

moving outward toward the periodic orbit. The cycles converge to the periodic orbit 
asymptotically as t ^ oo for VF^*, and r — >■ — oo for W^+. 

There are curves that leave move out to a maximum distance within the 
range ^/2 < w < \/3 {2M < r < 3M) when t; = 0, then turn around and move to p~ . 

We label these manifolds , and all of these geodesies have > L^. The union 
of these manifolds fill a region labeled Ti in Figured In the limit L ±|Lc|, the 
solutions asymptote to the union of and . The — > oo limit {h — > 0) are 

a rather interesting class of geodesies, which we discuss in the following subsection. 

Similarly, we have solutions that come in from negative infinity (3~) for w < 0, 
reach a minimum distance from the black hole in the range w > \/3 (r > 3M) when 
V = 0, then go back out to positive infinity (3^) with v > 0; they lie on invariant 
manifolds we call W^^, and are labeled T2 in Figure [21 As with Ti, all these geodesies 
have angular momentum squared > L^. In the limit L — >■ ±|Lc|, they asymptote 
to the union of the manifolds and W^± , the manifolds that connect 5~ to 
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and to 5+ respectively. The solutions on asymptotically (r — > oo) approach 

the circular orbit as they spiral towards it from and those on (beginning at 

T = —oo) asymptotically spiral away from the circular orbit to 3"*". 

The remaining two regions of the flow we call Ui and U2- Ui consists of the 
manifolds W^'^, which spiral in from 5^ into the black hole. They do not encounter 
a turning point in r. These geodesies all have < L^, and in the limit L — ^ ±|ic| 
asymptotically approach the union of manifolds W^oo and . Similarly, Ui consists 
of the manifolds W^, which spiral out from the white hole to 3^, do not have a 
turning point in r, have < L^, and in the limit L — > ±|ic| asymptote to the union 
of W^_^ and . 

The L'^ ^ oo (h^O) limit 

In region Ti, geodesies with h correspond to the limit = L^/E'^ — )■ oo. Fixing 
the physical energy E to be finite, these geodesies have several curious propertied, 
most easily deduced from the geodesic equations in Schwarzschild coordinates ([5]). 
First, for r < 2M, 

dt r-V2 

— = ±—= (40) 
dr ^r{2M - rY + P{2M - r) 

where the +(— ) sign corresponds to the part of the trajectory in the white (black) hole 
where f > (r < 0). Taking the limit £ — )• cjo, one gets that these geodesies projected 
onto the (r, t) plane in the Penrose diagram correspond to t — constant lines (and 
recall that inside the horizon t is a spacelike coordinate); i.e. they emanate from 
the white hole singularity, turn around at the intersection of the event and Cauchy 
horizons at r = 2Af , then continue to the black hole singularity. Next, we calculate 
how much cycling motion in <f> they execute. From the geodesic equations 

'^^± ' (41) 

dr ^r^ +rP{2M -r) 

again where +(— ) corresponds to the white (black) hole regions. Taking the limit 
£ — )■ oo, and integrating from r = to 2M in the white hole and back from r = 2M to 
in the black hole gives |A$| = 27r for the journey — these geodesies circle in azimuth 
exactly once going from the white to black hole singularity. Finally, we compute the 
total affine time Act = J da from 

da r3/2 

^ = ± - , ■ (42) 

dr E^r^ +P{2M -r) 

For finite E, in the limit £ — ?> oo, A(7 = 0! This rather bizarre result could be remedied 
by an infinite rescaling of cr — > (t(£o/^), where is some finite constant with dimension 
of length to make the scaling dimensionless. It is unclear what the physical significance 
such a rescaling is. 

Topology of the flow 

The solutions in the flow Ti that spiral out from the white hole, turn around at 
\/2 < w < and return to the black hole, trace out surfaces that are topologically 

§ in contrast to the £^ — >■ oo limit region of T2, which simply correspond to geodesies that pass by 
the black hole with an infinite impact parameter 
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equivalent to two cones smoothly joined together at u = 0, with vertices anchored in 
the critical points at ; see Figure [31 The flows Ui coming from 5~ into the 

black hole topologically form an (infinite length) open cone, with the open end at 
and the vertex at p^; similarly for the flows in U2- The geodesies in T2 that come in 
from 5~ and return to 3+ lie on families of surfaces that are topologically equivalent 
to (infinite) cylinders. 




Figure 3. Topology of the flow traced out by solutions within each of the 
four invariant manifolds Ti,T2,Ui and U2 (see Figure [2}. Note that the sample 
trajectories are schematic only and for visual aid; only those solutions that flow 
close to the critical points trace out multiple cycles, and these are clustered 
near w = \/3 (r = 3M). 

Extension of Solutions Through Collision 

Let x(^) be any solution of ([T8|) which ends in collision at ^ = 0. That is, x(^) — > 
as ^ — >■ 0+, where ^ — )■ 0+ means that ^ approaches where ^ > 0. It is proven in 
[11] that x(^) is branch regularizable if and only if 7 € 5, where we define the set S 
as follows: Let m,n be relatively prime integers, then S — {m/n\n odd, n > m > 0}. 
Recall that 7 = (1 + (3)~^,a = 2/3. In our case, a = 3 and 7 = 2/5 belongs to the set 

We define branch regularizable as follows: A solution x — W(^) of (US]) which either 
begins or ends in collision at ^ = ^* is branch regularizable at ^* if it has a unique 
branch extension at A branch extension is defined by considering two solutions 
Wi(^) and W2(C) of dUl), where Wi(0 ends in colhsion at C and W2(0 begins 
in collision at Then, W2(C) is a branch extension of Wi(^) if W2(C) is a real 
analytic continuation of Wi(^) at = ^* for 1^ in a neighborhood of ^ = £,* ■ 
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The fact that x(^) is branch regularizable imphes that there is a way that solutions 
can smoothly and uniquely be extended through p+ in backwards time and through 
p~ in forwards time. For example, let Z(^) = (x(^),y(^)) be a trajectory in phase 
space defined for ^ < that collides with at ^ = 0. Then there exists a unique 
extension of Z(^) for ^ > 0. This extension is real analytic as a function of ^ in a 
neighborhood of ^ = 0, and corresponds to a smooth bounce of P in {xi,X2) space. 
In the spacetime picture, collision with corresponds to the geodesic encountering 
the black hole singularity. The field equations of general relativity do not describe 
how spacetime can be extended beyond a singularity, and it is usually thought that 
a theory of quantum gravity is required to "resolve" the singularity. Nevertheless, 
one way to map the branch regularized extension to geodesic motion "through" the 
singularity would be to identify the black hole singularity with a white hole singularity 
of a second Schwarzschild solution of identical masf[i]. 

6. Conclusions 

In this paper we have studied the relationship between the null geodesic structure of 
the Schwarzschild black hole solution, and the corresponding inverse-cubic Newtonian 
central force problem, using the methods of McGehee. Both these problems have 
been well studied before, though what we believe is novel in this paper is highlighting 
the exact correspondence between the two descriptions, allowing insights from the 
dynamical systems approach to be brought to the geodesic problem, and vice-versa. 
Indeed, in that regard it is rather amusing to note that McGehee titled his paper 
"Double Collisions for Classical Particle System with Nongravitational Interactions" . 
It is also rather interesting that what in the Newtonian picture may be regarded as 
a clever "trick" using coordinate transformations to blow up the singular point of 
collision between two particles, is in a sense the natural way to describe Schwarzschild 
spacetime. Another remark is that understanding the invariant manifolds and unstable 
hyperbolic points allows standard techniques to be used to show that perturbations of 
the geodesic flow will generically cause chaotic motion. For example, for this purpose, 
solutions in Ti may be regarded as forming a homoclinic loop if we identify p^ and 
p~ . Perturbations will then generically break the homoclinic loop and cause chaotic 
motion by the Smale-Birkhoff theorem TU' . 

Even though we restricted attention to null geodesies for simplicity, we expect that 
similar mappings could be used for timelike particles, or more complicated geometries 
like Kerr. This would be an interesting avenue for future work. 
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